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ABSTRACT 



We study the process of spinning up black holes by accretion from slim disks in a wide range of accretion rates. We 
show that for super-Eddington accretion rates and low values of the viscosity parameter a (< 0.01) the limiting value of 
the dimensionless spin parameter a* can reach values higher than a, = 0.9978 inferred by Thorne (1974) in his seminal 
study. For AI — lOMEdd and a = 0.01 spin equilibrium is reached at a, = 0.9994. We show that the equilibrium spin 
value depends strongly on the assumed value of a. We also prove that for high accretion rates the impact of captured 
radiation on spin evolution is negligible. 

Key words, black holes physics — accretion disks 



1. INTRODUCTION 

Astrophysical black holes (BHs) are very simple objects - 
they can be described by just two parameters: mass and 
angular momentum. In isolation, BHs conserve the birth 
values of these parameters but often, e.g. in close binaries 
or in active nuclei, they are surrounded by accretion disks 
and their mass and angular momentum change. Accretion 
of matter always increases BH's irreducible mass and may 
change its angular momentum, usually described by the 
dimensionless spin parameter a* = a/M = J/M'^. The sign 
of this change and its value depend on the (relative) sign 
of accreted angular momentum and the balance between 
accretion of matter and various processes extracting BH's 
rotational energy and angular momentum. 

The question about the maximal possible spin of an ob- 
ject represented by the Kerr solution of the Einstein equa- 
tion is of fundamental and practical (observational) inter- 
est. First, a spin > 1 corresponds to a naked singular- 
ity and not to a black hole. According to the Penrose cos- 
mic censorship conjecture, naked singularities cannot form 
through actual physical processes, i.e. singularities in the 
Universe (except for the initial one in the B i g Ban g) are 
always surrounded by event horizons (jWaldl . Il984f ). This 
hypothesis has yet to be proven. 

In any case , the "third law" of BH thermodynamics 
ijBardeen et al.L [iQT^) asserts that a BH cannot be spun- 
up in a finite time to the extreme spin value a* = 1. 
Determining the maximum value of BH spin is also of prac- 



tical interest because the efhciency of accretion through 
a disk depends on the BH's spin value. For example, for 
the "canonical" value a* = 0.998 (see below) it is about 
77 fti 32%, while for g.^ 1 one has rj 42%. Recently, 
iBanados et al.l (|2009[ ) showed that the energy Eqom of the 
center-of-mass collision of two particles colliding arbitrary 
close to the BH horizon, grows to infinity (i?coM oo) 
when a* — 1^. 

A definitive study of the BH spin evolution will only be 
possible when reliable, non-stationary models of accretion 
disks and jet emission mechanisms are established. For now, 
one has to use simplified analytical or numerical models. 

iThornel (jl974l ) used the model o f a radiatively effi- 
cient , geometrically thin accretion disk (jNovikov &: Thorna , 
I1973D to evaluate BH spin evolution taking into account 
the decelerating impact of disk-emitted photons. The max- 
imum value so obtained, a* = 0.9978, has been regarded as 
the canonical value for the maximal BH spin. In this work 
we generalize Thome's approach, using models of advec- 
tive, optically thick accretion disks ("slim disks") to calcu- 
late maximum BH spin values for a large range of accretion 
rates. We show that for sufficiently large accretion rates 
they differ from the canonical value. 



^ Of more fundamental interest is the fact that the proper 
geodesic distances D between the ISCO and several other spe- 
cial Keplerian orbits relevant t o accretion disk struc ture tend to 
infinity _D — ^ oo when a* — >■ 1 (|Bardeen et al.l . Il972l ). 
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We start with a short discussion of previous work de- 
voted to the BH spin evolution. In Sect. [2] we give formu- 
lae for a general tetrad of an observer comoving with the 
accreting gas along the arbitrary photosphere surface. In 
Sect. [3] we give basic equations describing the BH spin evo- 
lution. The following section describes the applied model of 
slim accretion disks. In Sect. [S]we present and discuss the 
terminal spin values for all the models we consider. Finally, 
in Sect. [S] we summarize our results. 



1.1. Previous studies 

A number of authors have studie d the evolution o f the BH 
spin resulting from disk accretion. iBardeenI (jigTOT ) initiated 
this field of research by stating the problem and solving 
equations describing the BH spin evolution for accretion 
from the marginally stable orbit. Neglecting the effects 
of radiation he proved that such a process may spin-up 
the BH up to a* = 1. Once the classical models of ac- 
cretion disks were formu lated (jShakura fc Sunvaevl . Il973t 
iNovikov fc Thorn^ . 119731 ) , it was possible to account prop- 
erly for the decelerating impact of radiation (frame drag- 
ging makes counter-rotating photons m ore likely to be cap- 
tured by the BH). As mentioned above. [Thornel (|l974[) per- 
formed this study and obtained the terminal spin value for 
an isotropically emitting thin disk of a* = 0.9978, indepen- 
dently of the accretion rate. The original study by Thorne 
was followed by many papers, some of which are briefly 
mentioned below. 

T he first to challenge the universality of Thome's limit 
were |Ab ramowicz fc Lasotal fl980. 'l who showed that geo- 
metrically thick accretion disks may spin up BHs to termi- 
nal spin values much closer to unity than the presumably 
canonical a* = 0.9978. Their simple argument was based 
on models bv iKozlowski et al] (|1978D who showed that for 
high accretion rates the inner edge of a disk may be located 
inside the marginally stable orbit; with increasing accre- 
tion rate, arbitrarily close to the marginally bound orbit. 
However, this conclusion assumed implicitly a low viscosity 
parameter a, whereas for high viscosities the s ituation is 
more complicated (see lAbramowicz et all . 120101 and refer- 
ences therein). 

iModerski et al.l ([1998') assessed the impact of possible 
interaction between the disk magnetic field and the BH 
through the Blandford-Znajek process. They showed that 
the terminal spin value may be decreased to any, arbitrarily 
small value, if only the disk magnetic field is strong enough. 
Given the current lack of knowledge about the strength of 
magnetic fields (or the magnetic transport of angular mo- 
mentum in the disk, see iGhosh fc Abramowicd . 119971 ) and 
processes leading to jet emission, a more detailed study 
cannot be performed. The situation may further be compli- 
cated by energy extracti on from the inner parts of accretion 
disks (jLivio et al.l . ll999[ ). 

iPopham fc Gammid (|1998D studied the spinning-up of 
BHs by optically thin advection dominated accretion flows 
(ADAFs). They neglected the contribution of radiation to 
BH spin as such disks are radiatively inefficient. They found 
that the terminal value of BH spin is very sensitive to the 
assumed value of the viscosity parameter a and may vary 
between 0.8 and 1.0. iGammie et al.l (|2004l ). besides mak- 
ing a comprehensive summary of different ways of spinning 
up supermassive BHs, presented results based on a single 



run of a GRMHD simulation (with no radiation included) 
obtaining terminal spin a* = 0.93. 

Cosmological evolution of spins of supermassive BHs 
due to hierarchical mergers and thin-disk accretion 
episodes has b e en re cently intensively studied. Although 
IVolonteri et al.l ()2005[ ) arrived to the conclusion that accre- 
tion tend to spin-up BHs close to unity, as opposed to merg- 
ers which, on the average, do not i nfluence the spin evolu - 
tion, the fo l lowin g studies by e.g. Volonteri et al.l ()2007[ ): 
iKing et al.l (|200i); iBerti fc Volonteril (|2008^ showed that 
the situation is more complex, the final spin values depend- 
ing on the details of the hi story of the accretion events (see 
als olFanidakis et all . 1201 1[ ) . 

iLi et al.l ( 2005f ) included the returning radiation into 
the thin-disk model of Novikov fc Thorne and calculated 
the spin-up limit for the BH assuming the radiation cross- 
ing the equatorial plane inside the marginally stable orbit 
to be advected onto the BH. Their result (a* — 0.9983) 
slightly differs from Thome's result, thus showing that re- 
turning radiation has only a slight impact on the process of 
spinning-up BHs. In our study we use advective, optically 
thick solutions of accretion disks and account for photons 
captured by the BH in detail. However, we neglect the im- 
pact of the radiation returning to the disk on its structure. 



2. THE TETRAD 

We base this work on slim accretion disks, which are not 
razor-thin and have angular momentum profile that is not 
Keplerian (for details on the assumptions made and the 
disk appearance see Sect. Therefore, photons are not 
emitted from matter in Keplerian orbits in the equatorial 
plane and the classical e xpressions for photon momenta 
(e.g., iMisner et al.l . I1973D cannot be applied. Instead, to 
properly describe the momentum components of emitted 
photons, we need a tetrad for the comoving observer instan- 
taneously located at the disk photosphere. Below we give 
the explicit expression for the components of such a tetrad 
assuming time and axis symmetries. A detailed derivation 
is given in Appendix [XI 

Let us choose the following comoving tetrad. 



(1) 



where is the four- velocity of matter (living in [t, (f>, r,9]), 

Nl is a unit vector orthogonal to the photosphere ([r, 9]), 

unit vector orthogonal to {[t, </)]), 

is a unit vector orthogonal to u\ TV' and Kq {[t, cj), r, 9]). 
The tetrad components are given by. 



ar 



-1/2 



_^ 9rr 

gee 



d6>, 
dr 



d6>, 
dr 



n -1/2 



-1/2 



u = 



rf + + vSl 



y/gtt + ^94,4,1^ - 2a;) 



(2) 



(3) 



(4) 
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^{l + A^v^)-^/^{Avu' + Sl), (6) 

where 6 = 9^.{r) defines the location of the photosphere, rji 
and are the Killing vectors, I = u^/ut, — u'^/u*, oj is 
the frequency of frame-dragging and the expressions for v 
and SI are given in Eqs. (jA.lOp and (|A.4p . respectively. 

3. SPIN EVOLUTION 

3.1. Basic equations 

The equations describing the evolution of BH dimensionless 
spin parameter a* with respec t to the BH en ergy M and 
the accreted rest-mass Mo are (iThornd . [TqtI . 



dJ/M2 1 M0U4, + (^) 



rad 



dh^M = dlnM = MM^ut + {^)^^^ 



2a* 



dM 
dM) 



Ut 



(AM 



/Mo. 



(7) 



(8) 



rad 



The energy and angular momentum of BH increases due to 
the capture of photons according to the following formulae. 



(rfA/)rad= / T'^VkNdS, 
J disk 

(dJ)rad - / r^'^kN^dS, 
J disk 



(9) 



(10) 



where rjk and ^fc are the Killing vectors connected with 
time and axial symmetries, respectively, T'^^ is the stress- 
energy tensor of photons, which is non-zero only for photons 
crossing the BH horizon and dS*, the "volume element" in 
the hypersurface orthogonal to , is given by Eq. IB. 
From Eqs. ^ and ^ it follows that 



/dM 
I "dT 



27r /Tout 



rad -J 
dj\ r^- /•'■out 



rad 







where 



dS = d(j)dr [g^^ - guQ^^, 



1/2 / /d^* 

grr +909 [-^ 



(11) 



(12) 



(13) 



3.2. Stress energy tensor in the comoving frame 
Let us choose the tetrad given in Eq. ([T]): 

elo^^u' el)=N^ (14) 

pi — pi — qi 

^{2) — '5(3) — 

The disk properties, e.g., the emitted flux, are usually given 
in the comoving frame defined by Eq. ([T^ . The stress ten- 
sor components in the two frames (Boyer-Lindquist and 
comoving) are related in the following way, 



=T(")('3)e|„)ef^). (15) 
The stress tensor in the comoving frame is 

T(^m^2 / loSCn^^^Ti'^f^Uinhdhdb, (16) 

where loS — Io{r)S{ii, b) is the intensity of the emitted ra- 
diation, d and b are the angles between the emission vector 
and the and S^ vectors, respectively, C is the capture 
function defined in Sect. 13. 4[ the factor 2 comes from the 
fact that the disk emission comes from both sides of the 
disk and tt^") = are the normalized components 

of the photon four-momentum in the comoving fra me. The 
latter are given by the following simple relations ijThorneL 

^(0) = 1, 
TT^^' = COS a. 



(17) 



TT*-^-* = sin a COS 6, 



r(3) 



sin d sin b. 



Eqs. (HI]) and (HI]) take the form. 



/dM 
I dt 



rad J disk 



T 



e\^)elp)ilkNAS, (18) 



^) = I T^-'^^'^e^.^fUWS, (19) 

/ rad "'disk 



where e'^^^ is our local frame tetrad given by Eq. p4 
Taking the following relations into account, 



we have, 



(20) 



(a) 



(21) 



l'{p)Tlk 



where. 



(24) 
(25) 



Therefore, Eqs. (fT5)) and ([T^ may be finally expressed as, 
/dM\ _ /■^'^ r^'"^ 



47r 



\ 'l^ / rad "'r-i„ -'0 

X cosasinadadby^dr. 



A J\ /•'•out /•27r ^7r/2 



J Tin Jo Jo 

X COS d sin d dd db y/^ dr, 



with 



fT _ f 2 \ 1/2 / , f d()* 

'9 = [9t<i> ~ 9tt g4,4> ) \ 9rr + gee I 



(26) 



(27) 



(28) 
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3.3. Emission 

The intensity of local radiation may be identified with the 
flux emerging from the disk surface, 

lo = F{r). (29) 

The angular emission factor S is given by (lThornel[T97l . 

f /tt isotropic 



S{a,b) I (3/77r)(l + 2 cos a) limb darkening ^"^^^ 

for isotropic and limb-darkened cases, respectively. In this 
work we assume that the radiation is emitted isotropically. 



3.4. Capture function 

The BH energy and angular momentum are affected only b y 
photons crossing the BH horizon. Following iThornj (|1974l) . 
we define the capture function C, 



C 



I if the photon hits the BH, 
in the opposite case. 



(31) 



He rein, we c alcula te C in two ways. First, we use the origi- 
nal iThornd ()1974l) algorithm modified to account for emis- 
sion out of the equatorial plane. For this purpose we calcu- 
late the constants of motion, j and k, for a geodesic orbit 
of a photon in the following way. 



J = a. 



(32) 



k = 



[ttI - {n^ + aMTTtsmO^f I SIT? e^] (33) 



which replaces Thome's Eqs. (AlO). This approach does 
not take into account possible returning radiation, i.e., a 
photon hitting the disk surface is assumed to continue its 
motion. Such a treatment is not appropriate for optically 
thick disks - returning photons are most likely absorbed or 
advected towards the BH. 

To assess the importance of this inconsistency we adopt 
two additional algorithms for calculating C. Using photon 
equations of motion we determine if the photon hits the 
disk surface ((Bursa, 2006). Then we make two assumptions, 
either the angular momentum and energy of all "return- 
ing" photons are advected onto the BH (Ci), or all are 
re-emitted carrying away their original angular momentum 
and energy, and never hit the BH (C2). In this way we es- 
tablish two limiting cases allowing us to assess the impact 
of the returning radiation. 

We note here that for fully consistent treatment of the 
returning radiation (as in iLi et al.l . l2005l for geometrically 
thin disks) it is not enough to modify the capture function, 
solving for the whole structure of a self-irradiated accre- 
tion disk is necessary. The latter has not yet been done for 
luminous and geometrically thick disks. We are currently 
working on implementing such a scheme and will study its 
impact onto the BH spin evolution in a forthcoming paper. 



lAbramowicz et al l ('1996') and 'Ga mmie fc PophamI (|1998D . 
Here, we follow [Kato et al. (2008) and assume the poly- 
tropic equation of state when performing vertical integra- 
tion. The forma l ism w e use in this section was adopted from 
ISadowski et all (120111 ). 

In the structure equations we take G = c = 1 , and make 
use of the following expressions involving the BH spin: 



A = r2 



2Mr 



A = + + 2Mra^ , 



C = I- 3r-' 
V ^ 1- 2r-^ 

n = 1- 4a*r; 



o -3/2 

2a*r* 
3afrl 



(34) 



with a* = a/M and 
We also define 



r/M. 



(35) 



and dimensionless accretion rate m — M / M-^dd, where 
-^^Edd — 16iEdd/c^ is the critical accretion rate approxi- 
mately corresponding to the Eddington luminosity (iEdd = 
1.25 X 1O^^M/M0 erg/s) for a disk around a non-rotating 
black hole. 

The equations describing slim disks, written in the cylin- 
drical coordinates are: 
(i) Mass conservation: 



il = -27rSAi/2. 



1,^ 



(36) 



where S ~ p dz is disk surface density, while v denotes 
gas velocity as measured by an observer co-rotating with 
the fluid and is related to the four-velocity u'" by Ru^ — 

(ii) Radial momentum conservation: 



dv 



where 



A 



A 1 dP 

1-v^dR ^ T.dR' 

MA {n - - 
R^An+n- 



1 - n^R^ 



(37) 



(38) 



and il — u'l' /u* is the angular velocity with respect to a sta- 
tionary observer, J7 = fi — w is the angular velocity with re- 
spect to an inertial observer, Q.^ = ±M^^'^ j (R?^'^ ±aM'^^'^) 
are the angular frequencies of the co-rotating and counter- 
rotating Keplerian orbits and R — Aj (R? A}-^"^) is the radius 

of gyration. P = J^j^ pdz is the vertically integrated total 
pressure. 

(iii) Angular momentum conservation: 



M ^ ^ Ai/2Ai/2r 



(39) 



where C = u^,, £{„ is a cons tant and F is the Lorentz factor 
(jGammie fc Pophaml . [19981 ) : 



4. SLIM ACCRETION DISKS 

4.1. Equations 

In this section we pr esent sl i m dis k equations. They were 
derived originally bv lLasot"al (|l994f ) and improved e.g., by 



F2 = 



(iv) Vertical equilibrium: 



2o2 



A 



{2N + 3) 



P 



(40) 



(41) 



4 



A. S^dowski et al.: Spinning up black holes with super-critical accretion flows 



(v) Energy conservation: 



M 



P d In P , 



E dlni? 



PdlnS 



P d In 773 
S dlni? 



d In 7^4 
dlni? 



where the amount of heat advected Q^'" is: 



Iff- 

Assuming the polytropic index N 



- 3 we have: 



"Til 
V2 

m 



^ [ T^dz^ 128/315 iJ 

-|r / pTdz = 40/45 
^Jo Jo 



(42) 



(43) 



(44) 



fc40 



P V5/3- 1^45 
- [ pz^dz = 1/18 



256 
315' 



The equations given above form a two-dimensional system 
of ordinary differential equations with a critical (i.e., sonic) 
point. For each set of disk parameters the regular solution 
is possible only for one specific value of which is an 
eigenvalue of the problem. The appropriate value may be 
found using either the relaxation or the s hooting method . 
For details on the numerica l procedures see lSq,dowskil (|2009[ ) 
and lS£|.dowski et al.l (|201lh . 

4.2. Disk appearance 

In this section we briefly describe the properties o f slim disk 
solutions. For a more detailed discus sion see e.g.. ISadowskil 
(|2009[ ). lAbramowicz et af] ()2010[ ). or lBursa et all (|201lD . 

The radial profiles of the emitted flux for a non-rotating 
BH are presented in Fig.[TJ For low accretion rates (to ^ 1) 
they almost coincide with the Novikov & Thorne solutions 
(the small departure is due to the angular momentum taken 
away by photons, an effect which is neglected in our slim 
disk scheme) . When the accretion rate becomes high advec- 
tive cooling starts to play a significant role and the emission 
departs from that of the radiatively efficient solution. This 
departure is visible as early as for to = 1 — the emission 
extends significantly inside the marginally stable orbit. For 
super-critical accretion rates the flux increases monoton- 
ically towards the BH horizon. Different colors in Fig. [T] 
denote solutions for different values of the viscosity param- 
eter a. Although they are very similar, one can notice that 
the higher the value of a, the lower the accretion rate at 
which advection starts to modify the emission profile. 

In Fig. [2] we plot disk thickness profiles (cos Qh = H/r) 
for a range of accretion rates and two values of a. For 
to > 0.1 the inner region of the accretion disk is puffed 
up by the radiation pressure and the disk surface corre- 
sponds to the location where the radiation pressure force 
(proportional to the local flux of emitted radiation) is bal- 
anced by the vertical component of the gravity force. For 
the Eddington luminosity (to ss 1) the highest H/R ratio 



10^' 



' — ' 10^' 

\ 

e 

o 102^ 
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fn = 100.0 


a- 0.01 : 






a= 0.1 ^ 


: m = 1 


oTo~'^S^- 




r / / 
1 / 


fn = liT^^^^^; 






/m = 0.1~~~\^ 




\ \ 1 

/ / 


1^ = o.oT~"^-\^^ 




- 1 \ 
' / ' 


\ 




\ 


. ISCO 






10 


100 




r [M] 





Fig. 1. Flux profiles for M = 10 M© and a* = 0.0. 



0.001 




Fig. 2. Photosphere profiles profiles for M — 10 M© and 
a* = 0.0. 



equals ^ 0.3 (cos 0// « 0.3), while for the largest accretion 
rate considered (to = 100) it reaches ^ 1.5 (cos0/f « 0.83). 

In the thin disk approximation the angular momentum 
of gas follows the Keplerian profile. This condition is not 
satisfied for advective accretion disks with significant ra- 
dial pressure gradients. In Fig. |2] we present angular mo- 
mentum profiles for disks with different accretion rates, 
CL = 0.01 (left) and a = 0.1 (right panel). It is clear that 
the higher the accretion rate, the larger the departure from 
the Keplerian profile. However, the quantitative behavior 
depends strongly on a. For a < 0.01 the fiow is super- 
Keplerian in the inner part (e.g., between r = 4.5A/ and 
r — 1AM for m = 100). For larger viscosities (a > 0.1) 
and high accretion rates the flow is sub-Keplerian at all 
radii. As a result, the value of the angular momentum at 
the BH horizon (£in) also depends strongly on a, decreas- 
ing with increasing a. Dependence of the flow topology on 
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_m=1.0 



m=0.01, 0.1 



a= 0.01 

Keplerian 




12 14 16 18 20 




10 12 14 16 IB 20 



r [M] 



r [M] 



Fig. 3. Profiles of the disk angular momentum for a 
The spin of the BH a* = 0. 



0.01 (left) and a = 0.1 (right panel) for different accretion rates. 



th e viscosity paramet e r has been recently studied in detail 
by lAbramowicz et al.l (|2010h . 



5. RESULTS FOR BH SPIN EVOLUTION 

Using the slim disk solutions described in the previous 
section we solve Eqs. (O and ^ using regular Runge- 
Kutta method of the 4th order. To calculate the inte- 
grals (Eqs. pip and (112ll )jwe use the alternative extended 
Simpson's rule (jPressl . 1200211 basing on 100 grid points in 
a, b and radius r. We have made convergence tests proving 
this number is sufficient. 

In Figs. |3]and[S]we present the BH spin evolution for 
a = 0.01 and 0.1, respectively. The red lines show the 
results for different accreti o n rat es while the black line 
shows the classical 'Thorn^ (|1974 ) solution based on the 
[Novikov & Thorne (1973) model of thin accretion disk. Our 
low accretion rate limit does not perfectly agree with the 
black line as the slim disk model does not account for the 
angular momentum carried away by radiation. As a re- 
sult, the low- luminosity slim disk solutions slightly over- 
estimate (by no more than few percent) the emitted flux 
lea ding to s t ronge r deceleration of the BH by radiation — 
the iThornel (|1974[) result is the proper limit for the low- 
est accretion rates. When the accretion rate is high enough 
(e.g., 771 > 0.1), the impact of the omitted angular momen- 
tum flux is overwhelmed by the modification of the disk 
structure introduced by advection. 

It is clear the spin evolution is quantitatively differ- 
ent for different values of the viscosity parameter. For the 
lower value (a = 0.01) the BH spin can reach values sig- 
nificantly higher than 0.998 for the highest accretion rates 
(a* = 0.9995 for m = 100), while for the higher viscosity 
{a = 0.1) the terminal spin value decreases with increasing 
accretion rate down to a* — 0.9800 for m = 100. This be- 
havior is connected with the viscosity impact on the critical 
point topology. Generally speaking, the higher the value of 
a, the lower the angular momentum of the fiow at BH hori- 



0.9999 




m=100 (a. = 0.9995) 

m=10 (a. = 0.9994) 

m=l (a. = 0.9988) 

m = 0.1 (a. = 0.9967) 

m = 0.01 (a. = 0.9966) 

Thorne (a. = 0.9978) 



0.00 



Fig. 4. Spin evolution for a = 0.01. 

zon for a given accretion rate (Fig. [3]) , leading to a slower 
acceleration of the BH rotation. 

To study this fact in detail we calculated the rate of BH 
spin-up for "pure" accretion of matter (without accounting 
for the impact of radiation). For that case the BH spin 
evolution is given by (compare Eq. (O): 



da^ _ 1 
dlnAf ~ M ut ~ 



(45) 



In Figs.[n]and[7]we plot with black lines the first term on the 
right hand side of the above equation for different accretion 
rates and values of a. The red lines on these plots show the 
absolute value of the second term. The intersections of the 
black and red lines denote the equilibrium states i.e., the 
limiting values of BH spin for pure accretion. These val- 
ues differ significantly from the previously discussed results 
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Fig. 5. Spin evolution for a = 0.1. 



only for low accretion rates. In contrast, for high accre- 
tion rates radiation has little impact on the spin evolution 
and the value of terminal spin is mostly determined by the 
properties of the flow. In Fig. [5] we plot the radiation im- 
pact parameter ^, defined as the ratio of the disk-driven 
terms on the right hand sides of Eqs. ([7]) and (|45L 



(46) 



If the captured radiation significantly decelerates the BH 
spin-up, this ratio drops below unity. On the other hand, it 
is close to unity for BH spin evolution which is not affected 
by the radiation. According to Fig. [8] the latter is in fact 
the case for the highest accretion rates independently of a. 

In Table [T] we list the resulting values of the terminal 
BH spin for all the models considered. The first column 
gives results for our fiducial model (A) including Thome's 
capture function, emission from the photosphere at the ap- 
propriate radial velocity. 

The second column presents results obtained assuming 
the same (Thome's) capture function and profiles of emis- 
sion, angular momentum and radial velocity as in the model 
A, but assuming the emission takes place from the equa- 
torial plane instead of the photosphere. The resulting ter- 
minal spin values are equal, up to 4 decimal digits, to the 
values obtained with the fiducial model. This result is as 
expected for the lowest accretion rates, where the photo- 
sphere is located very close to the equatorial plane. For the 
highest accretion rates the location of the emission has no 
impact on the BH spin-up, as the spin evolution is driven 
by the fiow itself and the effects of radiation are negligi- 
ble. However, for moderate accretion rates one could have 
expected significant change of the terminal spin. We find 
that the location of the photosphere has little impact on 
the resulting BH spin regardless of the accretion rate. 

Our third model {V) neglects the flow radial velocity 
when the radiative terms are evaluated. Similar arguments 
to those given in the previous paragraph apply. For the 
lowest accretion rates, the radial velocity is negligible and 
therefore it should have no impact on the resulting spin. 



Fig. 6. The rate of spin-up or spin-down by "pure" accretion 
(radiation neglected) for a — 0.01. Profiles for five accre- 
tion rates are presented. Their intersections with the red 
line (marked with blue crosses) correspond to equilibrium 
states. For the two lowest accretion rates the equilibrium 
state is never reached (a, — > 1). 




Fig. 7. Same as Fig.[n]but for a = 0.1. 



For the highest accretion rates, the spin-up process depends 
only on the properties of the fiow. Once again, however, the 
impact of this assumption on moderate accretion rates is 
not obvious. The radial velocity turns out to be of little 
importance for the calculation of the terminal spin (only 
for a — 0.1 and moderate accretion rates the difference 
between models A and V is higher than 0.01%). 

In the fourth and fifth columns of Table [1] results for 
models with the same assumptions as the fiducial model, 
but with different capture functions are presented. The first 
alternative capture function (Ci) assumes that the angular 
momentum and energy of all photons returning to the disk 
are added to that of the BH. This assumption has a strong 
impact on the spin evolution — the terminal spin values are 
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fn = 100.0- 




Fig. 8. Radiation impact factor ^ (Eq. [IB]) for different ac- 
cretion rates and values of a. The dotted line corresponds 
to the thin-disk induced spin evolution. For ^ 1 spin evo- 
lution is not affected by radiation. Stars denote the equi- 
librium states (compare Table [T]). 



higher, sometimes approaching a* = 1. This may seem sur- 
prising because in the classical approach the captured pho- 
tons are responsible for decelerating the spin-up. This effect 
comes from the fact that the cross-section (with respect to 
the BH) of photons going "against" the frame dragging is 
larger than of photons following the BH sense of rotation. 
As frame dragging is involved, this effect is significant only 
in the vicinity of the BH horizon. For our model Ci, how- 
ever, the probability of photons returning to the disk does 
not appreciably differ for co- and counter-rotating photons, 
as they both hit the disk surface mostly at large radii. 

The other capture function (C2) assumes, on the con- 
trary, that all returning photons are re-emitted from the 
disk with their original angular momentum and energy (and 
never fall onto the BH). This assumption cuts off the pho- 
tons which would hit the BH in the fiducial model after 
crossing the disk surface. Thus, we may expect decreased 
radiative deceleration and increased values of the terminal 
spin parameter. However, these changes are not significant, 
reflecting the fact that most of the photons hit the BH di- 
rectly, along slightly curved trajectories. Only for a — 0.1 
and moderate accretion rates do the terminal spin values 
differ in the 4**^ decimal digit. 

Neither of the models with a modified capture functions 
is self-consistent. To account properly for the returning ra- 
diation one has to modify the disk equations by introducing 
appropriate terms for the outgoing and incoming fluxes of 
angular momentum and additional radiative heating. No 
such model for advective, optically thick accretion disk has 
been constructed. The emission profile should be signifi- 
cantly affected (especially inside the marginally stable or- 
bit) by the returning radiation leading to different rates of 
deceleration by photons. In view of o ur result s for m odels 
C\ and C2, as well as the results of iLi et all (|2005[ ). one 
may expect the final spin values for super-critical accretion 
flows to be slightly higher than the ones obtained in this 
work. 



The last column of Table[T]gives terminal spin values for 
"pure" accretion (radiation neglected) . Under such assump- 
tions the BH spin could reach a* = 1 for sub-Eddington 
accretion rates as there are no photons which could de- 
celerate and stop the spin- up process. As discussed above, 
for the highest accretion rates the resulting BH spin val- 
ues agree with the values obtained for the fiducial model as 
radiation has little impact on spin evolution in this regime. 

6. SUMMARY 

We have studied BH spin evolution due to disk accretion 
assuming that the angular momentum and energy carried 
by the fiow and the emitted photons is the only process af- 
fe cting the BH r otation. We generalized the original study 
of lThorn^ (jl974 ) to high accretion rates by applying a rela- 
tivistic, advective, optically thick slim accretion disk model. 
Assuming isotropic emission (no limb darkening) we have 
shown that 

(i) the terminal value of BH spin depends on the accre- 
tion rate for to > 1, 

(ii) the terminal spin value is very sensitive to the as- 
sumed value of the viscosity parameter a — for a < 0.01 
the BH is spun up to a, > 0.9978 for high accretion rates, 
while for a > 0.1 to a* < 0.9978, 

(iii) with a low value of a and high accretion rates, the 
BH may be spun up to spins significantly higher than the 
canonical value a* — 0.9978 (e.g., to a, — 0.9994 for a — 
0.01 and to — 10) but, under reasonable assumptions, BH 
cannot be spun up arbitrarily close to a» — 1, 

(iv) BH spin evolution is hardly affected by the emitted 
radiation for high (to > 10) accretion rates (the terminal 
spin value is determined by the flow properties only), 

(v) for all accretion rates, neither the photosphere pro- 
file nor the profile of radial velocity significantly affects the 
spin evolution. 

We point out that the inner edge of an accretion disk 
cannot be unique ly defined for a super-critical accretion 
(jAbramowicz et al.. , . 2010.1 , as opposed to geometrically thin 
disks where the inner edge is uniquely located at the 
marginally stable orbit {Rms)- In the thin-disk case the 
BH spin evolution is determined by the flow properties at 
this particular radius (as there is no torque between the 
marginally stable orbit and BH horizon) and the profile of 
emission (terminating at Rms)- For super-critical accretion 
rates, however, one cannot distinguish a particular inner 
edge which is relevant to studying BH spin evolution. On 
the one hand, the values of the specific energy (ut) and the 
angular momentum (u^) remain constant within the stress 
inner edge. On the other, the radiation is emitted outside 
the radiation inner edge. These inner edges do not coincide 
as they are related to different physical processes. 

Our study was based on a semi-analytical, hydrodynam- 
ical model of an accretion disk which makes a number of 
simplifying assumptions like stationarity, no returning radi- 
ation, aP prescription and no wind outflows. Thus, the re- 
sults obtained in this work should be considered only qual- 
itative, as the precise values of the terminal spin parameter 
are very sensitive (e.g., through viscosity) to the flow and 
emission properties. Moreover, we neglected the impact of 
magnetic fields and jet ejection mechanisms. Nevertheless, 
our study shows that Thome's canonical value for BH spin 
(a, = 0.9978) may be exceeded under certain conditions. 
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Table 1. BH spin terminal values 



capture function 




C 


C 


C 


Ci 


C2 




model: 






A 


T 


V 


A 


A 


NR 


thin disk 






0.9978 


0.9978 


0.9978 


0.9981 


0.9978 


1 




rh = 


0.01 


0.9966 


0.9966 


0.9966 


0.9975 


0.9966 


1 




m — 


0.1 


0.9967 


0.9967 


0.9967 


1 


0.9967 


1 


a = 0.01 


rh - 


= 1 


0.9988 


0.9988 


0.9988 


1 


0.9988 


0.9998 




rh = 


10 


0.9994 


0.9994 


0.9994 


-> 1 


0.9994 


0.9996 




rh = 


100 


0.9995 


0.9995 


0.9995 


1 


0.9995 


0.9995 




rn = 


0.01 


0.9966 


0.9966 


0.9966 


0.9975 


0.9966 


^ 1 




rh — 


0.1 


0.9975 


0.9975 


0.9975 


1 


0.9975 


^ 1 


a = 0.1 


rh - 


= 1 


0.9924 


0.9924 


0.9923 


^ 1 


0.9927 


0.9948 




rh = 


10 


0.9846 


0.9846 


0.9845 


0.9901 


0.9847 


0.9951 




rh = 


100 


0.9800 


0.9800 


0.9800 


0.9803 


0.9800 


0.9801 



Notes. C - Thome's capture function, Ci - all returning photons adverted onto the BH, C2 - all returning photons neglected; A 
- our fiducial model, T - emission from the equatorial plane, V - zero radial velocity, NR - pure accretion, radiation neglected. 
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Appendix A: The tetrad for an observer 

instantaneously located at the photosphere 

Our aim is to derive the tetrad of an observer moving along the photo- 
sphere that would depend only on the quantities which are calculated 
in accretion disk models, i.e., on the radial and azimuthal velocities 
of gas and the location of disk photosphere. 

The metric considered here is the Kerr geometry in the Boyer- 

Lindquist coordinates [t,(f>,r,9]. The signatur e adopted is H . 

Similarly as in Carter's Les Houches lectures llCarterl . 119721 ). we will 
consider two fundamental planes; the symmetry plane So = [t, 0] and 
the meridional plane A^* = [r,0]. (Four)- Vectors that belong to the 
plane So, will be denoted by the subscript 0, and vectors that belong 
to the plane Al*, will be denoted by the subscript *. For example, 
the two Killings vectors are jJq, ^q. Note, that for any pair Jfg, one 
has. 



A.l. Stationary and axially symmetric photosphere 
A. 1.1. The photosphere 

Numerical solutions of slim accretion disks provide the location of 
the photosphere given by -ffph('') = rcos9. This may be put into 
r cos 8 — Hpij(r) = F{r, 6) = 0. The normal vector to the photosphere 
surface has the following [r, 6] components. 



Nl = Af* 



dF dP 



NL 



de, 

— , ] 

dr 



where 



de,{r) 
dr 



dF 

dr , 



IdF 



(A.l) 



(A.2) 



is the derivative of the angle defining the location of the photosphere 
at a given radial coordinate [cos0«(r) = H-p\i{r) /r\. Its non-zero com- 
ponents after normalization [(A'^*A', ) = — 1] are 



Nl 



de* 

dr ^~ 



-1/2 



Qrr 



de. 

"dr" 



-1/2 



(A.3) 



9 



A. S^dowski et al.: Spinning up black 



holes with super-critical accretion flows 



n; = i-gee) 



-1/2 



1 + 



Qrr ( dS, 

dr 



-1/2 



There are two unique vectors S, confined in the [r, 9] plane which are 
orthogonal to Af« (and therefore are tangent to the surface). From 
{S,Nt) = and (St, St) = — 1 one obtains the non-zero components 
of one of them: 



1 

Qrr 



1 fde. 

Quo V dr 



-1/2 



(A.4) 



Ad, 
dr 



1 

Qrr 



1 



aeo V dr 



AO, 



-1/2 



A. 1.2. The four-velocity of matter and the tetrad 

The four- velocity u of gas moving along the photosphere may be de- 
composed into 



where 



"0 = Ao {v' + ^e) 



(A.5) 



{A.6) 



is the four-velocity of an observer with azimuthal motion only. The 
normalization constant Aq comes from {uquo) = 1 and equals 



Aq = [gtt + ^g4>4.{p.-2u)\ 



{A.7) 



It is useful to construct a spacelike vector (kq) confined in the [i, <f>\ 
plane, that is perpendicular both to u and uq. From (kk) = — 1 and, 
e.g., {kuq) = we have. 



{A.8) 



where I = u^/ut is the specific angular momentum. Note that the set 
of vectors [uq, 7V^, Kq, SX\ already forms the desired tetrad valid for 
the pure rotation (n"" = 0) case. 
The normalization condition (uu) = 1 gives, 



A = [gtt + f^S0</,(r2 - 2w) - v 



,21-1/2 



{A.9) 



where v is related to the radial component of the gas four- velocity 
by: 

2 K/gn' (gtt + ng^^jQ - 2w)) 

V = = . (A.IU) 

1 + (u-/S:f 

The vectors we have just calculated (n, kq) are both orthogonal to Nt 
since {NtSt) = 0. To complete the tetrad we need one more spacelike 
vector (S) that is orthogonal to these three. Let us decompose it into. 



5^ = au^ + I3k'o + -fNl + SSi 



(A.ll) 



The orthogonality conditions (koS) = and (NtS) = give immedi- 
ately 7 = /9 = 0. The only non-trivial condition is (uS) = 0. Together 
with (SS) = -1 it leads to: 



= (1 + A2t,2)-l/2(^^u» + S^J^ 



{A.12) 



A. 2.1. The four velocity 



Similarly as in Section IXOl we may always uniquely decompose 
a general timelike unit vector, into 



= A{u^ + vSl), 



(A.14) 



where is a timelike unit vector, and SI is a spacelike unit vector. 
Formula HA. 141 uniquely defines the two vectors Uq , SI and the two 
scalars A, v. The vectors and scalars 



{A,v,ul„Sl}, 



(A.15) 



can be calculated from known quantities given by slim disk model 
solutions. 

The four- velocity IIA.141 1 defines also the instantaneous 3-space of the 
comoving observer with the metric ■yn^ and the projection tensor hj^, 



lik = gik -UiU^., 

h\ = 5\ - u' Uk 



(A.16) 
(A.17) 

We define two unit vectors Kq and Nl by the unique condition, 

(ko«o) = 0, {S,N,) = 0. (A.18) 

As before, the four vectors 

^\a) = K. 4: K, Si], (A.19) 

which form an orthonormal tetrad of an observer with the four- 
velocity Uq can calculated from the solutions of the slim-disk equa- 
tions. 

A. 2. 2. The photosphere 

In the most general case of a non-stationary and non-axially symmet- 
ric photosphere, the location of the photosphere may be described by 
the following condition, 

F{t,(f>,r,e) = 0. (A. 20) 

The vector N normal to the photosphere has the components. 



dF dF dF 8F 



(A.21) 



which may be calculated from slim disk solutions. 
Let us project N into the instantaneous 3-space of the comoving ob- 
server 1 IA.I7I 1 and normalize to a unit vector after the projection. 



|(AfAf)|i/2' 



(A.22) 



In terms of the tetrad I IA.191 1. such constructed vector has the 
following decomposition. 



TV* = TV [a K) 1 (a:) -h 7 (4) + -5 (5:)]. 

The components N,a,^,& are known. 



(A.23) 



spacetime: 



= [u\ Nl, Kl, S']. 



{A.13) 



This tetrad is known directly from the slim disk solutions, as it de- 
pends on the calculated quantities {u^, O, I and 6*{r)) only. Any 
spacetime vector X^, could be uniquely decomposed into this tetrad 
with, XfA) = Xi e»j^^. 



A. 2. The general case 

In this section we will assume nothing about the four- velocity of mat- 
ter and the location of photosphere. Both may be non-stationary 
and non-axially symmetric. Following the same framework as in the 
previous sub-Section, we will describe how to obtain the tetrad of 
an observer instantaneously located at the photosphere that depends 
only on the quantities calculated by accretion disk models. 



A. 2. 3. The tetrad 

Let us now write decompositions of the four vectors: the first two we 
have derived, the next two guessed (but the guess should be obvious): 



= A[l{ul) +0{NI) +0{k'o) +V{Sl)], 
= N[a{ul) + l{Nl) +-y{Kl) +5{Sl)], 
= k [0(4) +b{Nl) +1{kI)) +0{SI)], 
S' = S [A{ul) + B(Nl) + C{K,l)) + l{Sl)]. 



(A.24) 
(A.25) 
(A.26) 
(A.27) 



The four unknown components, b,A,B,C one calculates from the 
four non-trivial orthogonality conditions {{uu) = by construction, 
cf. 1 IX24I1 and iEjSJ), 



(uS) = 0, (NS) = 0, (5k) = 0, (Nk) = 0, 



(A. 28) 
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and the two unknown factors , 
ditions 

(kk) = ■ 



S, from the two normahzation con- 



-1, (SS): 



(A. 29) 



The condit ions IIA.2811, iKldl are given by Unear equations. 
Equations IIA.24t - llA.29ll define the tetrad e of an observer comov- 
ing with matter, and instantaneously located at the photosphere: 

e = [u', N', K.', S']. (A.30) 

Both the matter and the photosphere move in a general manner. The 
zenithal direction in the local observer's sky is given by N^. 



and 



Inserting this into the formula for dA we get dA 
gtt 94,4,)^^'^ ■ The final formula for dS is, 



(B.7) 



dS = dtd<Adr 



9?0 ■ 



- 9tt g4,4, 



1/2 



Qrr 



1 de, 

Vdr 



(B.8) 



Appendix B: Integration over the world-tube of the 
photosphere 

For stationary and axially symetric models we have so far, 

A' = Nl = unit vector orthogonal to the photosphere. It is in the 

[r, &\ plane 

S% = unit vector orthogonal to A' that lives in the [r, &\ plane 
ji' = four-velocity of matter. It lives in the [t, 0, r, 6] space-time 
= Kq = unit vector orthogonal to [/* that lives in the [t, </>] plane 

S" = unit vector orthogonal to , and ft' . It lives in the [t,4>,r,9] 
space-time 

e = [m% A', K*, S"] = the tetrad comoving with an observer 
located in the photosphere 

The integration of a vector over the 3-D hypersurface H or- 

thogonal to A'^* (i.e. the 3-D world-tube of the photosphere) may be 
symbolically written as 

/ i...)iN*dS, (B.l) 
Jn 

where dS is the "volume element" in H. 

Obviously, the hypersurface H is spanned by the three vectors 
[ji*, K*, S*]jv Each of them is a linear combination of 5', S^Jat, 
and each of the three vectors from [r?',^*, SI]m is orthogonal to A'. 

Therefore, one may say that the hypersurface H is spanned by 
5% S'JJjv- It will be convenient to write 



dS = dAdR, dR = dr^ grr + gee (-J^^ ' (^-2) 

where dR is the line element along the vector SI, i.e. along the pho- 
tosphere in the [r, 9] plane, with 6 = dt {t) defining the location of the 
photosphere, and dA is the surface element on the [t, </>] plane. 

In order to calculate dA, imagine an infinitesimal parallelogram with 
sides that are located along the t = const and <j> = const lines. The 
proper lengths of the sides are du = \gtt\^^^dt and dv = \g4,4,\^^^d(j> 
respectively, and therefore dA, which is just the area of the parallel- 
ogram, is given by 

dA = du dv sin a = dtd<f> \gtt\^^'^ \g4>4)\^^'^ sin a (B-3) 

where a is the angle between the two sides. Obviously, the cosine of 
this angle is given by the scalar product of the two unit vectors rii 
and Xi pointing in the [t, 4>\ plane into t and <^ directions respectively. 
These vectors are given by (note that rii = ZAMO), 



(V,t) (V,0) 



|3J'=(V,t)(Vfet)[V2' |9^HV,</')(Vfe0)|i/2- 
Because (V^t) = 5*^ and (Vi</>) = 5'^^, one may write. 



5*, 5 



(B.4) 



(B.5) 



|g"|l/2' \g4>4>\i-/2' 

Therefore, 
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3.8 



3.6 



3.4 



3.2 



1 1 1 1 1 1 

A 


II 1 II 


1 1 1 1 1 

A 


1 1 1 1 1 1 1 1 

a= 0.01 


A 




□ 


a= 0.1 


- (Thorne) 




---ED-- 
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III 1 1 II 


1 1 1 1 1 1 


□ 

1 II 
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